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FILM BOILING OF MERCURY DROPLETS 

by Kenneth J. Baumeister*, Glen J, Schoessow* 
and Charles E. Chmielewski* 


Vaporization times of mercury droplets in 
Leidenfrost film boiling on a flat horizontal plate 
are measured in an air atmosphere. Extreme care was 
used to prevent large amplitude droplet vibrations 
and surface wetting; therefore, these data can be 
compared to film boiling theory. For these data, 
diffusion from the upper surface of the drop is a 
dominant mode of mass transfer from the drop. A 
closed-form analytical film boiling theory is devel- 
oped to account for the diffusive evaporation. 
Reasonable agreement between data and theory is 
seen. 

Introduction 

Up to the present time, in the study of film 
boiling of liquid droplets, most experiments have 
considered fluids such as water, hydrocarbons, or 
cryogens. References 1 to 6 provide a comprehen- 
sive summary of past work plus a discussion of the 
phenomenon of droplet film boiling, often called 
Leidenfrost boiling. Recent interest in liquid- 
metal fast -breeder reactors and concern for their 
safe operation has lead to the need for film boiling 
deta for liquid metals. In this paper, the authors* 
report vaporization times of mercury droplets in 
film boiling on a flat horizontal plate in an air 
atmosphere and correlate the data. The heat trans- 
fer coefficient, plate temperatures, etc, are also 
reported. 

Foppendiek, et al, ^ave reported some film 
boiling data of mercury droplets. 7 *® They reported 
that large amplitude droplet vibrations and inter- 
action of the droplet with the surface occurred dur- 
ing the vaporization process. The authors report 
the vaporization time of only a single size droplet 
and unfort' jmtely give only an average droplet rad- 
ius rathe:- than an exact volume measurement. 

: eliminate the problems described by 

Poppendiek, et su.. only very small mercury droplet 
volumes (0.3? cm in diam. or less) were used herein 
to obtain stable film boiling with a minimum of 
droplet vibration and surface wetting. 

In an air atmosphere in the volume range con- 
sidered, (Diffusion from the upper surface of the 
drop is a dominant mode of mass transfer as shown 
by an evaluation of the dimensionless groups de- 
veloped herein ana shown by some experimental data 
from Ref. C. The previously-developed closed form 
theories for the heat transfer coefficient and va- 
poriz- tic- tines either neglected or considered 


diffusion from the upper surface of the drop as a 
small part of the mass transfer process. As a re- 
sult, these theories, such as developed in Refs. 1 
or 10, can not be used to correlate the data. The 
numerical scheme presented in Ref. 2 could be ap- 
plied. However, since a closed-form analytical 
theory is convenient, & closed- form solution was 
sought and is presented herein for the heat trans- 
fer coefficient and vaporization time of droplets 
in diffusion -dominated film boiling. 

The analysis is limited to two-component sys- 
tems involving a pure single -component liquid with 
small solubilities for the gas and a single- 
component gas of high purity (with the exception of 
air for which measured values of the diffusion co- 
efficient exist). 

Experimental Apparatus and Procedure 

The experimental equipment consisted basically 
of a heated flat surface insulated on its sides. 

The initial droplet vaporization experiments on the 
hot surface were conducted ir. a cylindrical tank 
60 cm in diameter and 60 cm long with & helium at- 
mosphere. This apparatus was the same one used in 
Ref. 8 to measure the effects of tne environmental 
diffusivity on the vaporization times of the drops. 
The first data point taken for a drop in a dry en- 
vironment was always repeated at the end of a data 
set to check for repeatability. 

Although the remote handling was inherently a 
safer way to conduct the vaporization experiments 
with mercury which is toxic, remote handling made 
close observation and accurate volume control dif- 
ficult. Instead of using the closed chamber, the 
flat heated plate was operated in an open air at- 
mosphere inside a hood. In this way, the mercury 
droplet could be placed by hand on the surface in a 
more gentle manner and thereby eliminate much of 
the droplet vibration problems discussed in Ref. 7 
and seen in our earlier experiments. 

The oxidatioii of the mercury drop did not ap- 
pear to be much of a problem based on the volume of 
the residue left at the end of the evaporation. In 
addition, a check was made against data in the hel- 
ium atmosphere (approximately the same theoretical 
diffusivity) and found to be in agreement. Also, as 
will be shown later, in either air or an inert at- 
mosphere, the dominant mechanism of mass transfer 
is diffusion. Air is the only gas for which a 
measured value of the diffusion constant for mer- 
cury is given. 12 
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Mercury vapor contamination of the room air was 
of concern even with the hot plate inside the ven- 
tilated laboratory hood. The initial check using 
an air flow meter indicated that instabilities of 
air under the ventilation hood in the forms of 
swirls and eddys could cause leakage of mercury vap- 
or into the room. A detailed ventilation study was 
performed. It was found that the experiment could 
be conducted safely with the hood window opened to 
25 cm and an air flow of 1 m/sec. All the pre- 
scribed government safety recommendations 1 * were 
followed including the taking of blood and hair 
samples. 


be considered in a later section. 

Previous Correlation of Vaporization Times 

Previous equations for predicting the vapori- 
zation times of drops over predicted the time re- 
quired for the mercury drops to vaporize by as much 
as 300 percent. This was surprising since these 
theories had correlated water, hydrocarbons and cry- 
ogens over a very large range of droplet volumes. 
Reference 10 had even correlated water vaporizing 
into a he" 1 ' »*m atmosphere where diffusion had a ma- 
jor contribution. 


The mercury was delivered to the heated sur- 
face using a thin pipette attached to a hypodermic 
syringe. The mercury was preheated for some of the 
experiments by an electrical resistance type heater 
wound around a metal storage cup. For these hot 
liquid measurements, the delivery pipette was 
stored with its tip in the hot mercury to reduce 
cooling during delivery. 

An electrical timer which measured time to one 
tenth of a second was used to measure the vaporiza- 
tion time. The vaporization times ran between 15 
and 60 seconds. The plate temperature was measured 
using imbedded chromel-alumel thermocouples and ex- 
trapolating to find the surface temperature. The 
volume of each mercury drop expelled onto the hot 
plate was found by measuring the length of the liq- 
uid column it had occupied in a pipette. The pi- 
pettes were calibrated using laboratory scales to 
develop calibration curves in terms of the length 
of mercury in the pipette (measured to the top of 
the miniscus). Prior to delivery, the pipette was 
placed next to a scale graduated in millimeters 
with a readability of 0.5 mm or an accuracy of 
±0. 00025 ml in volume. 

Data were J ^st taken with cold mercury, and 
then the mercury was heated to 265° C in a stainless 
steel cup. It was noticed that when volumes above 
0.014 ml (0.32 cm diem.) were used, the droplet 
would jump and vibrate violently when first deliv- 
ered onto the heated plate and this erruption of 
the droplet would usually occur again during evap- 
oration. Some of these droplets never reached a 
stable film boiling condition but continued to 
bounce across the plate and vibrate throughout the 
entire evaporation process. Consequently, drops 
with diameters less than 0.32 cm (0.014 cc) were 
used in the experiment. 


Diffusive evaporation from the upper surface 
of the aroplet explains why earlier correlations 
did not work on the mercury data. As shown in Ap- 
pendix A of this paper, if the following dimension- 
less group is greater than 2, diffusion from the 
upper surface of the drop will be a dominant mech- 
anism. 


Ndc 
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These and other symbols are defined in the Nomen- 
clature list at tM end of this paper. Since 
is greater than 2 for mercury drops, the previously- 
developed closed-form theories require some modifi- 
cation. 

Poppendiek^ vaporized mercury drops in film 
boiling into a saturated atmosphere of mercury va- 
por. His results are shown in Fjg. 2. When dif- 
fusion is minimized by reducing the partial pressure 
driving force, the vaporization times increase by a 
factor of 7. Clearly, both the theory and experi- 
ment indicate that diffusion from the upper surface 
of the drop can be dominant mechanism in mercury 
film boiling. 

Unfortunately, Poppendiek did not document the 
volume of mercury used so that the saturated at- 
mospheric data could not be compared to theory. 

In add! cion, as will be discussed later, sur- 
face wetting beneath the drop (nucleate boiling) 
and convection augmented diffusion from the upper 
surface of the drop could also contribute to the 
difference between previous theories and the data. 

Analysis of Vaporization Times 


Experimental Results 

The measured vaporization time data are dis- 
played in Fig. 1 for mercury Initially at room tem- 
perature and preheated to 2650 C. As seen in 
Fig, 1, preheating the mercury had a very small ef- 
fect on the vaporization time. Since the rutio of 
the heat of vaporization to C*AX required to raise 
the drop to its boiling point ft small, these re- 
sults were expected. 

The overall average heat transfer coefficients 
to the mercury droplets can be estimated from these 
curves; however, it is first necessary to estimate 
the area through which the heat passes. This will 
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The experimentally measured vaporization time 
of a discrete liquid drop undergoing film boiling 
can be determined by a direct integration of an en- 
ergy balance on the drop; 

2L . ^(vlAgWAr (2) 

dt 

where the total heat transfer coefficient, h^, and 
the drop* s bottom area A fi are dependent on the 
volume V of liquid that exists at any time, t. 

In order to integrate eq, (2), h^ and A, 
must be related to the droplet volume, properties 
of the liquid and vapor, plate temperature, and the 
environmental conditions surrounding the drop. The 


$ 




2 





expression for the heat transfer coefficient is de- 
veloped by analyzing the model shown in Fig, 3. 

Heat transfer to the drop occurs by conduction 
across the small vapor gap and by radiation from 
plate to the drop. Mass is lost by diffusive evap- 
oration from the upper surface of the sphere and 
by film boiling off the bottom of the drop. The 
complete analysis for h* is presented in Appen- 
dix A. The relationship of A fi to droplet volume 
is given in Appendix B. 


First, eq. (2) will be non-dimensional i zed us- 
ing the same definitions in Ref, 10. Let 
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Eq. (2) can now be written as 
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Integrating eq. (8) gives: 
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Substituting the expression for eq. (A15) and 

a£, (rr 2 eq. (B2)), into eq. (9) and performing 
tne Integration gives 


t * . 6«0 2 v^-664 

where 



(10) 


( 11 ) 
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Equation (10) will be compared to the experi- 
mental data in a later section. 

Heat Transfer Coefficient 


The various heat transfer mechanisms combine to 
give a total heat transfer coefficient that can be 
found by rearranging eq. ( 2 ) into the form 


4 * 


.W dV 

^ dt 






( 12 ) 


The dV/dt term was found at selected volumes by 
measuring the slope of the vaporization curves in 
Fig, 1. The Ag was evaluated using eq. (B2) in 
Appendix B. 

The heat transfer coeff?cient data are dis- 
played in Fig. 4. As seen in Fig. 4, the heat 
transfer coefficients are very large compared to 
the normal values of 0.01 cal/ cm sec c (or 50 BTU/ 
hr ft 2 F) expected from film boiling theory. The 
analytical predictions for the heat transfer coef- 
ficient fall below the experimental values. How- 
ever, the analytical trend seems to follow the ex- 
perimental trend. Consequently, if desired, an 
empirical correction factor could be applied to the 
analytical expression to give agreement over a wide 
volume range. 

Because the emissivity of mercury is small, 
the heat is transferred to the drop primarily by 
conduction across the vapor gap. Therefore, the 
vapor gap beneath the drop can be estimated by the 
equation 


b 



(13) 


The gap thickness is estimated to be on the order 
of one ten-thousand of a centimeter. The hot plate 
surface was polished initially to a 5 micro -cm fin- 
ish. 


There are at least two possible reasons why 
the analytical heat transfer coefficient under- 
estimates the measured heat transfer coefficient. 
First, convection currents above the drop could 
enhance the diffusive mechanism from the upper 
surface of the drop. Second, because of the very 
small vapor gap thickness and the large surface 
tension of mercury, some direct contact between 
mercury and the surface could be occurring. It 
should be expected that analytical models, which 
assume an insulating vapor film and quiescent 
drops, will underestimate the measured heat trans- 
fer coefficient. Bredfiel<A 5 indicated that 
liquid-solid contact can occur even with large 
vapor gap thickness (0.001 cm) and low surface ten* 
•ion fluids. Harvey^(p. 42) indicates that a large 
percent change in the droplet diameter can occur 
due to direct contact. In addition, Marangonl 
forces could also play a role ii, enhancirg the heat 
transfer coefficient. 
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Correlation of the Vaporization Time 

The theory, eq. (10), is compared to experiment 
in Pig. 5. Reasonable agreement is seen especially 
at the smaller droplet volumes. The deviation be- 
tween experiment and theory could be due to the ex- 
tra heat that is transferred to the drop by direct 
contact or removed by convection augmented diffusion 
from the upper drop surface as was just discussed. 

It is expected that solid-liquid contact will 
decrease at the smaller droplet volumes. In this 
range, the surface tension forces will compact the 
mercury into a tight sphere and thereby reduce the 
pos; bility of liquid stir face contact. In this low 
volt ’je range, the vaporization times ere seen to be 
in improved agreement with the measured data. 

Conclusions 

Diffusive evaporation from the upper surface of 
the drop is a major mechanism of mass transfer for 
liquid mercury drops in Leidenfrost film boiling in 
an air atmosphere. Theoretical expressions are de- 
rived for the heat transfer coefficient and the va- 
porization times, which are in fair agreement with 
experimental data. 
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Nomenclature 

bottom area of drop 

area of drop in which diffusion occurs 

dimensionless area, defined by eq. (5) 

specific heat at c instant pressure of 
vapor 

diffusion coefficient 

coefficient of gravity 

gravitational constant 

heat transfer coefficient 

film boiling heat transfer coefficient 
when radiation and diffusion are assumed 
zero 

radiative heat transfer coefficient 
convection heat transfer coefficient 
total heat transfer coefficient to drop 
thermal conductivity of vapor 
characteristic length, eq. (6) 
molecular weight 

dimensionless diffusion parameter defined 
by •«. (A4) 

dimensionless diffusion parameter defined 
by *q. (11) 
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t p 

T s 

£T 

t 

t* 

V 

V* 

w(6) 

Z 
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X 


a 

o 
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saturation pressure 

gas constant 

radial corrdinate 

dimensionless radius, r/L 

drop radius 

contact radius 

temperature 

plate temperature 

saturation temperature 

temperature difference, (Tp - T fi ) 

time 

dimensionless time, eq. (3) 
droplet volume 

dimensionless volume, eq. (4) 
axial velocity at bottom of drop 
axial coordinate 
vapor gap thickness 

vapor gap thickness when diffusion ana radi- 
ation are not present 


emissivity of liquid mercury (c^ 
latent heat of vaporization 
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modified latent heat of vaporization defined 
in Ref. 10 as 


c as 
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viscosity of vapor 
density of liquid 
vapor density 
surface tension 

Appendix A 

Hc»t Transfer CoefflcUnt 


The overall energy balance on the droplet shown 
in Fig. 3 can be written as 
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On the right-hand side of eq. (Al), the first 
term represents heat removed from the drop by evap- 
oration from the lower surface, while the second 
term represents the latant heat requirements for the 
evaporation by diffusion from the upper surface. 

For the form of diffusive energy loss, the partial 
pressure of the liquid in the vapor environment is 
assumed small and an equivalent sphex^cal shape has 
been assumed. Heat transfer to the drop occurs by 
conduction across the small vapor gap 6, by radia- 
tion from the plate to the drop and by natural con- 
vection from the hot vapor and atmosphere surround- 
ing the drop. 

The momentum equation along with the appropri- 
ate boundary conditions was solved in Ref. 12 to 
give 

2 1 (p L ’ *v>» 3 (Xd (A2) 

3 ul 


w(B) 


where V and A- are dimensionless volume and 
area of the bottom of the drop respectively. They 
are defined by eq. (4) and (5) in the body of this 
report. 

Equation (A2) is now combined with the energy 
balance, eq, (Al), to determine the gap thickness of 
the vapor 5, 
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and as shown in Ref. 12 for pure film boiling, no 
radiation or diffusion 
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and for the heat transfer coefficient 

where the latent heat of vaporisation has been re- 
placed by the modified latent heat of vaporisation 


(A7) 


which accounts fur convect.on effects in the vapor 
gap (see list of symbols;. 

The dimensionless group is & measure of 

the importance of diffusion from the upper surface 
of the drop to conduction beneath the drop. 

Equation (A3) can be written as 
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The exact solution of eq. (AS) is presenter in 
Fig. 6. For small values of 
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The solution for 6/B 0 can be written as 
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For large values of N. 
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The solution for the heat transfer coefficient 
and vaporization times presented in Ref. 10 was 
based on eq. (A10), that is, diffusion is relativ- 
ely small compared to heat conduction across the 
vapor gap. For mercury drops in either air or an 
inert atmosphere, eq, (A10) is not valid because 
N dc is large (in the range of b to 10 for the mer- 
cury data). The balance of this appendix is con- 
cerned with the derivation of the heat transfer co- 
efficient using eq, (A12) where diffusion dominates. 

The total heat transfer coefficient can now be 
found by suktituting from eq. (A12) into the left- 
hand side of eq. (Al). 
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This form of the heat transfer coefficient could be 
obtained directly from eq. (Al) by simply neglect- 
ing the effect of vapor flow (-p^AwA^) beneath the 
drop. The analysis has indicated that this is a 
permissible assumption if - h^/hg is great- 

er than 2. 

Non-dimensionalizing eq. (A13) is tne nanner 
suggested by eq. (7) in the body of the report 
gives 
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The expression for hj, depends on volume 
through the terms r Q , A d , and A fi . The volume de- 
pendence of r Q , A d , and A fi are given respectively 

by eqs. (B4), and (B6). Substituting these values 
into eq. (A14) gives 

^ - v* OS? 


30 . 

5 BT. 


/kV(p. 

I AT 


plate is neglected. For V greater than 5. the 
drop is assumed to be represented by a flat cylin- 
der, such that only the upper area contribute? to 
diffusive evaporation. 
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Area of Drop 

The theory of Ref. 10 based the area of heat 
transfer on the maximum radius of the drop. How- 
ever, -or dimensionless drop volumes less than 5, 
that is 

v* = v _ < 5 

r °s c v 7 * <»> 

the contact area is considerably smaller than the 
area subtended by r , This is shown pictorially 
in Fig. 7. The curves were generated from the data 
of Kef. 14. The curves can be fitted by the follow- 
ing equations: 
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As an analytical simplification, the equation for 
the drop’s surfaco area in eq. (?6) is assumed to 
be represented by a complete sphere. That is, the 
area dim to the flattening of the sphere at the 


6 






HEATED PLATE 

Figure 3. - Schematic model for evaporation of flat bottom sphere. 



DIMENSIONLESS VOLUME, V* 

Figure 4. - Dependence of heat transfer coefficient on dimensionless 
volume tor mercury droplets. 
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